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Abstract 

We point out that when the decay of one electroweak scale super-WIMP state to another occurs 
at second order in a super- weak coupling constant, this can naturally lead to decay lifetimes that 
are much larger than the age of the Universe, and create observable consequences for the indirect 
detection of dark matter. We demonstrate this in a supersymmetric model with Dirac neutrinos, 
where the right-handed scalar neutrinos are the lightest and next-to-lightest supersymmetric part- 
ners. We show that this model produces a super-WIMP decay rate scaling as m^/(weak scale) 3 , 
and may significantly enhance the fraction of energetic electrons and positrons over anti-protons 
in the decay products. Such a signature is consistent with the observations recently reported by 
the PAMELA experiment. 
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1. INTRODUCTION 



The detection of non-gravitational cold dark matter (DM) interactions is an important 
goal that drives some aspects of modern particle physics, cosmology and astrophysics. The 
well-measured abundance of dark matter [l| does not specify the origin or specific nature of 
DM. In particle physics, when the DM particle masses are comparable to the electroweak 
scale (m x ~ v) there are two distinct (well- motivated) possibilities. 1 The first one, is the 
weakly interacting massive particle (WIMP) framework, in which the initially thermalized 
abundance of WIMPs is reduced via their weak annihilation, starting at the particles freeze 
out temperatures Tj ~ 0.05 m x {2, 3]. The second possibility (the super- WIMP) postulates 
that the DM particles have an interaction rate much weaker than the weak interactions. In 
this case, one assumes that the initial cosmological abundance of such particles is small, and 
that the rate for their direct production from thermal SM states, controlled by the square 
of the super-WIMP coupling yswi remains small relative to the Hubble expansion rate 
throughout the entire history of the Universe. The ratio of the super-WIMP thermalization 
rate to the Hubble rate at the weak scale, T th /H ~ 10~ 4 yg W x (Mp;/weak scale), gives a 
crude estimate for the resulting number densities of elecroweak scale super- WIMPs weighted 
by entropy. Observations require this ratio be ^ 10~ 12 , which sets the benchmark value for 
the super-WIMP coupling, y sw < 10~ 12 . The stability of electroweak scale super- WIMPs 
would have to be ensured by some parity in the dark sector. 

While the WIMP framework gives hope to the goal of direct detection of DM 0, 2, [f| (see 
however [7]), electroweak scale super- WIMPs would necessarily have a very tiny interaction 
strength that generally cannot lead to direct detection. Conversely, indirect signatures of 
DM, such as energetic gamma-rays, positrons and anti-protons, may be created by either 
WIMPs or super- WIMPs, provided that the latter decay with a lifetime longer than the age 
of the Universe. 

The recent claim by the PAMELA collaboration 2 jjjj] of an excess in the positron fraction 



above 10 GeV, is broadly in agreement with expectations for WIMP annihilation 15 



16 



13, 



1 Of course, it is also possible that nature could correspond to the parameter space between two extremes 
of WIMPs and super- WIMPs and/or have multiple DM states of each type. 

2 As well as the more recent claim by the ATIC collaboration of an electron excess at higher energies 
and other past hints of energetic electron/positron excess in the galaxy at high energies [9 , 
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18l . ligL [20|, |2l( or WIMP / super- WIMP decay [22j, |23| leading to an excess postron flux above 
the background of secondary postron production. Even though an astrophyiscal origin of 
this signal (unrelated to dark matter) is not ruled out, there have been numerous attempts 



to 
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ink the PAMELA positron excess to dark matter. Various analyses of the positron flux 
3, 22! find that the PAMELA result can be fit to models of dark matter, provided 
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that the annihilation and/or decay rate to positrons satisfy the following criteria: 



WIMP annihilation 
WIMP/super-WIMP decay 



(av) e+ ~ 0(3 x KT^cmV 1 ) x ( 
r e+ ~ O(10~ 51 GeV) x 



77K 



500 GeV 



1 TeV 



(1) 
(2) 



where m x is the mass of decaying/annihilating particles. One should keep in mind that these 
estimates have uncertainties both due to their dependence on the modeling of the propaga- 
tion of positrons in the galaxy (see 28|, |29| for a discussion) and due to uncertainties in the 
local dark matter energy density. In addition, PAMELA has also reported the measured 
flux of anti-protons 30J, which is well described by standard astrophysical production of 



anti-protons 



3l| 



While the annihilation rate in Eqn. ([I]) of WIMPs inside our galaxy is naively in conflict 
with the abundance-derived rate of (o~v) T f — 3 x 10~ 26 cm 3 s~ 1 , it has been argued that this 
may not necessarily be the case 32j, |33|, |34J , as significant enhancement factors are possible 



due to resonance annihilation and/or v~ l enhancement from the long-range attraction in 
the dark matter sector 3 . This can boost the galactic annihilation rate far above (av)^- As 
no excess in the anti-proton flux is reported, the positron excess suggests mostly leptonic 
channels of annihilation and small rates of decay to quark- ant iquark pairs for the DM. 

As a "boost factor" of size 10 — 10 3 needed for the interpretation of the PAMELA results 
via WIMP annihilation may be accomplished in many ways through model-building, an 
explanation via the decays of WIMPs or super- WIMPs may seem more ad hoc. In order 
to obtain the decay width of 10~ 51 GeV in Eqn. (EJ, one would typically introduce a decay 
constant x/decay ~ 10~ 25 — 10~ 23 , such that T decay oc y decay x (weak scale)/167r 2 . Such a 
coupling is twenty orders of magnitude smaller than a typical WIMP coupling and more than 
ten orders orders of magnitude smaller than the upper bound on the super- WIMP coupling. 
Although possible, such an option is not appealing, yet in most of the decay scenarios 
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discussed in the literature, such a coupling is introduced by hand to fit the PAMELA data. 

In this note, we consider the generic possibility when there are more than one super- 
WIMP states present in the DM sector, and the decay of one state to another is kinematically 
allowed. Such a decay can happen only in the second order of the super-weak coupling 
constant, and can be naturally suppressed down to the level given by Eqn. (j2j): 

(superWIMP)i -> (superWIMP) 2 + SM particles, T ~ 0{y A sw ). (3) 

The forth order of the superweak coupling, 0(yg w ), is necessary a tiny number not far from 
10 -50 , which may lead to a natural realization of the O(10 -51 GeV) decay width. 

We demonstrate this scenario realized in a supersymmetric Standard Model (MSSM) 
with exact i?-parity conservation extended with right-handed (RH) neutrino superfields and 
Dirac (or nearly Dirac) neutrino masses. In this model one identifies the Yukawa couplings 
in the neutrino sector y v ~ 10~( 12 " 13 ) with y$w 1^6 , 47, 48|. Analyzing the model, we find 
that the excess of energetic positrons and electrons may naturally arize from the decay of 
one RH sneutrino species into another z>^ — > v\£ + £~, and that for certain domains of the 
parameter space the production of anti-protons is inhibited. This particular realization of 
the decaying super- WIMP scenario is the subject of this paper. 

2. NATURALLY UNSTABLE SUSY DM CANDIDATES 

The spectrum of supersymmetric particles depends on the nature of supersymmetry 
breaking generating the soft masses of the MSSM Lagrangian. Currently, there is no over- 
whelming reason to adhere to any particular scheme of SUSY breaking; the only objective 
restriction on the SUSY-breaking mass pattern comes from the resulting phenomenology. 
Thus, the charged lightest supersymmetric particle (LSP) is ruled out due to the desire for 
a natural DM candidate, and a nearly exact i?-parity is required to keep the LSP stable if 
it is to be the DM. Many WIMP candidates such as left-handed sneutrinos are disfavored 



49j due to direct detection constraints. 

However, in order to accommodate the observed neutrino oscillations 4 the SM should 
be supplemented with right handed neutrinos v l R . The corresponding MSSM contains right- 



501 ] for a review of the experimental evidence of neutrino oscillation and [51( for a recent discussion 



See 

of the theoretical formulation of neutrino oscillation. 



4 



handed neutrino superfields that include both the fermionic right-handed neutrinos and their 
scalar partners, the right-handed sneutrinos v R . Further, any choice of scale for the Majorana 
mass for the RH neutrinos is technically natural, which leaves the Yukawa couplings in the 
neutrino sector a free parameter that can be varied in a wide range 10~ 13 < y l u < 1. If 
neutrino masses are Dirac (or nearly Dirac) then the Yukawa couplings are close to the 
lower end of this range, and the v R obtain the dominant contribution to their masses from 
the soft SUSY breaking mass terms in the MSSM; two of the v % R masses, resulting mostly 
from these mass terms, could be the LSP and the NLSP. 

Super- WIMP DM made up of v R was investigated recently in 



47 



481 ] and found to 



be a viable option. 5 Indeed, v R are in some ways a very promising candidate for DM in 



the galactic halo: they can have the right relic abundance [46|, |47|, |48|, and represents cold, 
neutral, colourless particles that do not interfere either with primordial nucleosynthesis or 
stellar evolution. Obviously, i> R DM is consistent with limits on the self-interactions of DM 
and leaves no detectable signal in direct detection searches 53|, |54j as their interactions are 
suppressed by y v . 

Despite these facts, one could question the phenomenological consequences of light Dirac- 
like neutrino masses due to the required smallness of the yukawa coupling. The SM Yukawa 
couplings already have a hierarchy of 10 6 between the smallest and largest Yukawas. This 
hierarchy, although puzzling, does not lead to problems with technical naturalness as chiral 
symmetry is restored in the limit of vanishing Yukawa couplings. The Yukawa couplings 
for Dirac neutrinos are very small O(10~ 13 ), which leads to an even larger Yukawa cou- 
pling hierarchy but are likewise technically natural. Further, Dirac neutrinos may not have 
lepton number violation which disfavors some forms of leptogenesis 6 as candidate theories 
for generating the observed matter-antimatter asymmetry of the Universe. However, many 
alternatives for baryogenesis exist such as electroweak MSSM baryogenesis 561], the Affleck- 



Dine mechanism [57] . or an effective modification of the EW phase transition [58l. l59l. 160 1 . 

The benefits of i> R as DM arguably outweighs the theoretical costs for the reason men- 
tioned in the introduction. If the PAMELA data is to be explained through the decay of a 
particle, then the decay width must be smaller than the width corresponding to the age of 



5 Decays of sneutrinos through R-parity violation was considered as a source of th e PAMELA excess in 52 1 . 

6 Note that leptogenesis with dirac neutrinos can occur in some models, see 55| for example. 
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the universe IV = r^ 1 ~ 10~ 42 GeV. Such small widths are challenging to naturally produce 
for weak scale WIMPs. On the other hand, in our scenario the small decay width is linked 
to another very small number, namely yl ~ (10 -13 ) 4 , which is an interesting possibility. 7 
One can immediately see that the decay widths are such that Y v <^iTu. These decays would 
still be occurring in the galactic halo leading to a primary source of positrons that could be 
the excess observed by PAMELA, and possibly explain the ATIC data as well. 



A. Dirac v l R as LSP and NLSP 



Following the initial stage of cosmological evolution, all superpartners of active SM 
species (charged under gaug e g roups) decay to v % R on the time scales controled by T oc 



y1 x (weak scale), see Refs. 47], |48] for details. For our scenario, it is sufficient to assume 
that the population of the NLSP (D R ) is not parametrically small compared to the LSP (i> R ) 
due to a very mild hierarchy of Yukawa couplings. 

Consider the i?-parity conserving MSSM supplemented with three Dirac v l R with the 
superpotential 

W = y*j e af3 H ua Lpp v R /i e a/3 H ua H df3 + ... (4) 

Here H ua = (H+,H®) and H da = {H%H^) are the up and down type Higgs chiral super- 
fields and Lp = (i>i,£2) is the lepton chiral superfield and e±2 = 1. We have not written the 
terms in the superpotential containing the quark and charged lepton Yukawa couplings as 
they are not directly relevant to our discussion. The Yukawa interactions are derived from 
the superpotential directly using 

pi 



Writing out the relevant component field terms in the resulting Lagrangian one has 

tyuk = v> Hi - P R P Hi) - n fa H~ - Hi flj) +h.c+-~ (6) 

When the lowest lying supersymmetric particles are (y l R , i> R ) the decay of one of these 
species into another can proceed via v R — > v R v k v l or v R — * v R l k i l . As we insist on exact 



7 This observation was made in [46| , where it was pointed out that the lifetime of the excited state of RH 
sneutrinos exceeds the lifetime of the Universe. 
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i?-parity, decays of v % R to purely SM particles is forbidden, while all other supersymmetric 
final states are kinematically inaccessible. 




FIG. 1: The tree level decay process for D R — > u^l^i^, v R — > i> R vv and the local operator 
approximation. 

The leptonic decays v R — > v R i k £ l proceed through Higgsino exchange 8 and the rate can 
be easily calculated as the function of m 1; m 2 and /i, while masses of the SM leptons can be 
safely neglected. The amplitude is given by 

iA = i{yyyT)U^i— 2 L\ (7) 

Here L is the left-handed lepton doublet, 1, 2 and j, I are flavour indicies of RH sneutrinos 
and SM leptons respectively, and p is the momentum carried by the Higgsino. Since all mass 
parameters are essentially free apart from the fi > mi > vri2 constraint, after of quoting a 
general formula, we present the answer in several different limits. If the Higgsino mass fi is 
larger than the sneutrino mass scale, the Higgsino propagator can be contracted, as shown 
in Figure 1, and the decay is mediated at leading order in p//i by the effective Lagrangian, 

£c ff = ^f- ffidpPti V^L\ (8) 

which is simply a product of sneutrino and left-handed lepton currents. Introducing effective 
couplings y\ = Ivl^ 2 , y\ = \uV\ 2 an d the energy release parameter A = m\ — m 2 

8 We disregard the issue of mixing between Higgsinos and gauginos as the mixing cannot introduce more 
than 0(1) corrections. 
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one can calculate the general decay width in a straightforward manner. One finds the result 
for the decay width given by 

T l+ A (24 m\ - 84Am? + 104A 2 mf - 50 A 3 mj - 8 A 4 m? + 20A 5 m 2 - 8A 6 mi + A 7 ) 
y\y\ 6144 m 3 7r 3 /i 4 



12 m 4 log (mi - A) 



vl 



(9) 

6144 m 3 vnV 

This general result obscures the physics somewhat, however, recall that in our scenario 
the sneutrinos are the LSP and the NLSP so necessarily the Higgsino is a larger mass scale. 
Let fi = nrrii where n is an order one number (n > 1) for this reason. Further, for this decay 
to occur m 2 < mi so let m 2 = xm± where x < 1, then expanding in small x one obtains 

r> = Z'ff 2 , (1 " 8 x 2 - 24 x 4 log (x) + 8x 6 + 0(x 7 )) . (10) 
ol44?7, 4 7r ci 

Thus we see that a good order of magnitude approximation for this decay width is IV/- ~ 
10~ 6 y\ y\ m\ for the whole range of parameter space. 

We also note that this decay width has a number of phenomenologically interesting limits. 
For example, considering A to be somewhat small on the scale of the sneutrino masses, we 
derive the total decay rate into charged leptons, 

This formula can be obtained directly from the neutron beta decay rate upon setting 
m e ,gA — ► 0, and identifying Gf cos9 c / y/2 — > 1/(4 /j 2 ). In another kinematic regime, when 
H is very close to m± and the energy release remains small, the rate is given by 

^ = 192^2 for A* - mi « A « mi, 2 . (12) 

Finally, if the scale of m 2 is small and can be neglected, while /x remains close to mi we have 
a third limit, 

= 1024 7T 3 for i" - ^' m 2 «mi. (13) 
We do not require that any of the hierarchies between mass parameters mi, m 2 , \i used to 
determine the widths r"+ 6 J- is strickly realized in our scenario. We have determined 
using these limits on mi, m 2 , \x to determine some phenomenologically interesting limits of 
the leptonic decay width. In what follows we will use the approximation T/+ j- ~ 10~ 6 yf y\m\ 
which is a good approximation due to the physics of our scenario, (ie a Higgsino heavier 
that the LSP and NLSP sneutrinos and m 2 < mi). 

8 



C. Hadronic v R decays as a source of p flux 

Depending on kinematics, besides the pure leptonic decay channels, there could be other 
decay channels potentially leading to energetic anti-quarks and eventually to anti-protons. 
Quark-antiquark pairs cannot be produced in the sneutrino decay via the intermediate gauge 
and/or Higgs bosons. The relevant part of the scalar potential that could lead to such decays 
is 

v = 2fi ytj f> R & H° d - 1{ h+) + £ vl k (yl k YH° u (# U °)X4 + h. c + ... (u) 

k 

These interactions induce three- and four-body decays of RH neutrinos: v R — > u R h(H), 
v R — > v R W + W~(H + H~), v R — * u R ZZ(AA), decays. To simplify our calculations, we shall 
assume that tan/5 = (H®) / (H®) is somewhat large, and the A, H, H ± Higgs bosons are 
heavy, so that only the decay to the lightest Higgs boson h and pairs of gauge bosons are 
kinematically possible. 

We concentrate on the three-body decay v R — > v R h, which is driven by the last term in 
Eqn. f[T4"l) when one of the Higgs is given a vev. The resulting decay width is 

2 12 v 2 ( . ml - A 2 ' 



r. = \y 2 | 2 x v h A + -± . (is) 

In this formula, Vh is the velocity of the outgoing Higgs boson, v = 246 GeV is the SM Higgs 
vev, and y\ 2 is the sneutrino-flavor changing combination of the Yukawa couplings: 

3 

If mi ~ v, Vh ~ 1 and y\ 2 ~ yiyi, the two-body decay rate would dominate over the 
leptonic rate by three orders of magnitude simply because of the phase space suppression 
of Eqs. fllip - ffTHl) relative to Eqn. (fl~5"j) is quite significant. The Higgs boson produced in the 
decay fragments further to gauge bosons and heavy quarks and leptons. For SUSY models 
its typical mass precludes it from decaying directly to W + W~ or ZZ, and bb decays are 
expected to dominate. Decays of this form is often encountered in neutralino annihilation, 
and the corresponding yield of antiprotons from a 6-quark injection has been evaluated 



61 



621 ] . The yield of antiprotons in the hadronization and subsequent decays of 66-pairs 
is quite significant, exceeding 10% per annihilation/decay event with a typical electroweak 
scale energy injection. Assuming an 0(0.1 — 1) yield of antiprotons in the v R v R h decay, 
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the ratio of Eqn. fll5l) to a typical three-body rate gives the relative strengths of antiprotons 
and positrons at the injection: 

^ ~ 10 3 x >< (17) 



A recent theoretical study of the antiproton fraction [31| limits the "boost factor" for the 

annihilation of dark matter to be less than 6(40) for the annihilation of 0.1(1) TeV particles. 

If we tune the model to fit PAMELA flux for positrons, this would translate into a constraint 

on Eqn. ffTTl) to be smaller than 0.1-1. Consequently, one would have to require a suppression 

of Vh\yf 2 \ 2 / {yf yV) down to the level of 10 -4 — 10~ 3 . Such a suppression may come from the 

smallness of the flavor-changing coupling yf 2 , or may originate kinematically, from small 

scale splitting, A < m/j. In both cases one would be able to fit PAMELA results without 

overproducing antiprotons. 

Note that decays that are loop suppressed will generally also have the difference in flavour 

decays to quarks and thus protons and anti-protons compared to leptonic decays. This is 

due to the exact R parity of our scenario and the fact that only the LSP v\ is kinematically 

accessible for the NLSP v x R to decay too. Loop decays with the same flavour structure do 

exist where the produced leptons form a loop that produces aZ/7 that subsequently decays 

to hadrons. The most problematic of these decays are suppressed by 

tf\ 22 
P ~ 10 3 x 9 %- x BR(Z -> hadrons) (18) 



$ e+ 167r 2 (4cos 2 (#w0) 
and the suppression is smaller than the required 0.1-1 



D. Numerology 

Our main result, EqB. (jllj) - (jl3j) can give an excellent fit to the PAMELA data in agreement 
with the requirement of Eqn.(j2]). Adopting T1+1- = 10~ 6 x y\y\ x m\ as an approximation 
for different kinematic regimes of the three-body decay rate, we conclude that the Yukawa 
couplings must satisfy the following relation: 

(2/ij/ 2 ) 2 ~10- 82 =>► (yiy 2 ) 1/2 ~ 1 x 10" 13 . (19) 

Are these values for the Yukawa couplings in agreement with Dirac neutrino Yukawa 
couplings suggested by the measurements of neutrino oscillations? The flavor basis in 
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the RH sneutrino sector and the active neutrino sector in general do not coincide and a 
one-to-one connection between neutrino phenomenology and Eqn. ffT9~]) is not possible. The 
flavor- changing combination of Yukawa couplings in Eqn. dlSp . remains uncertain and is not 
necessarily related to the measured mixing angles. If further assumptions are made about 
the structure of the PMNS and the Vr mass matricies in a particular constrained MSSM 
then the the flavour structure we have determined can be related to the PMNS mass ma- 
trix, but this is beyond the scope of this work. No immediate connection between our 
determined flavour structure and the obervable masses and mixing angles in the neutrino 
sector is possible without further experimental input, but we can check the consistency of 
Eqn. ([T9"j) with general expectations for the Yukawa couplings in the Dirac neutrino sector. 
Recall that Dirac neutrinos obtain their masses via m u = y u (H®) = y u v sin (3 / \/2 ~ y v x 175 
GeV, as we assume tan (3 to be large. To estimate the possible size of the Yukawa couplings 
we use the experimental determinations of the neutrino mass splitting. The difference in 



masses squared of the neutrinos have been measured by the KamLAND [63[ and K2K [64] 
collaborations to be 

[5m 2 ] atm ~ 2.8 x 10" 3 eV 2 , [5m 2 ] solar ~ 7.9 x 10" 5 eV 2 (20) 

and the absolute bound on the sum of neutrino masses from WMAP is given to be 

J2 m l = 0.67eVat95%CL. (21) 

which coincides with the limit on the heaviest neutrino mass m u < 0.7eV. 

In the case of the standard hierarchy of neutrino masses the heaviest neutrino is given 
by ~ (Sml) at m and one finds that the largest Yukawa coupling is 

/ m 2 \ 1/2 

y v ~ 3.0 x 10~ 13 v — j , (22) 

y V2.8 x 10~ 3 eV 2 / V ; 

which is perfectly consistent with the requirement on (y 1 y^) 1 ^ in our scenario. Alternatively, 
for the small neutrino mass splitting scenarios one could have Yukawa couplings close to the 
value suggested by possible given cosmological bounds, y v ~ 10~ 12 . In either case we see 
that the Yukawa couplings can realize our super- WIMP decay scenario that produces the 
PAMELA, and possibly ATIC, signals for electroweak scale DM masses. 
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3. DISCUSSION AND CONCLUSIONS 



We have presented an observation that the decay rates of super- WIMPs with the rate 
about ten orders of magnitude slower than the Hubble rate may originate from the transition 
of one super- WIMP state into another in the second order in the super- WIMP coupling so 
that T ~ y$ w x (weak scale). The resulting decay rates are just becoming possible to probe 
through indirect detection of dark matter via its decay products, such as electrons, positrons 
and antiprotons. 

Perhaps the most interesting aspect of our model is that super- WIMP physics is ex- 
tremely sensitive to the superpartner mass spectrum due to the existence of a very small 
coupling ysw- If 5 f° r example, one has a mass spectrum such that < m H i g gsino < m u 1 R ^ 
the sequential decays v\ — > Higgsino — > v\ will happen with r oc yl x (weak scale) oc 10 Hz, 
which is almost instantaneous on the scale of Eqn.([2]). A simple modification of the spec- 
trum to m^i < m^2 < mniggsino wa s shown to leads to an enormous delay in the v x R — > v\ 
decay. This delay results in an overall drop in the decay rate by 25 orders of magntiude, and 
makes the heavier component of the RH sneutrinos only very weakly unstable. The num- 
bers, given a huge disparities of different scales involved, work remarkably well, producing 
T e + ~ m^/(weak scale) 3 , which is in agreement with a putative explanation of PAMELA 
(and ATIC) signals by the decaying super- WIMPs. Moreover, for a moderate splitting be- 
tween two RH sneutrino components and/or an accidental suppression of the different flavor 
transition, the hadronic decays of v\ will be suppressed, which will in turn suppress the 
antiproton flux created by super- WIMPs. A kinematic suppression of the antiproton flux 
would imply a rather small energy splitting, A < rrih ~ 0(100) GeV between two super- 
WIMP states. Alternatively, if both PAMELA and ATIC are to be explained in the same 
way, one would have to choose A close to a TeV and assume a hierarchy \yf 2 \ 2 ^ v\v\ i n 
order to suppress the antiproton flux. 

One feature of our scenario may look somewhat unusual. We do not assume a complete 
degeneracy of RH sneutrino masses, and in fact requite them being split by the energy 
intervals comparable to their masses. This is not possible in the context of universal sfermion 
masses often employed in the SUSY literature, as any additional splitting induced by y v is 
minuscule. Therefore, one would have to imagine some additional theoretical mechanism to 
create such a splitting. 
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Although the example given in this note is quite natural, it is hardly unique. For example, 
one could envisage sequential i?-parity preserving decays of other super- WIMPs featured in 
the SUSY literature, such as gravitinos, axinos, modul(inos) etc. 

The decaying super- WIMP scenario as a tentative explanation of PAMELA and ATIC 
differs from the boosted WIMP annihilation in many ways. The most notable distinction is, 
of course, the early cosmology. The enhanced annihilation may lead to the nuclear-chemical 



consequences such as the overproduction of e Li [65(, and to extra ionization and a diffuse 
gamma background produced by the annihilating WIMPs. In contrast, the decaying 
DM scenario is immune to these potential problems. On the other hand, both annihilation 
and decay scenarios are subject to gamma ray and synchrotron emisson constraints (for a 



recent paper see 



671 ] ) coming from the central region of the galaxy where the density of dark 



matter is significantly enhanced. 
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